Stationary untrapped boundary conditions in 
general relativity 



00 ■ Roh-Suan Tung 

O ' 

' Center for Astrophysics, Shanghai Normal University, Shanghai 200234, China 

E-mail: tung@sh.nu . edu . cn 

rj^ ■ Abstract. A class of boundary conditions for canonical general relativity are 

, proposed and studied at the quasi-local level. It is shown that for untrapped or 

£NJ ' marginal surfaces, fixing the area element on the 2-surface (rather than the induced 

2-metric) and the angular momentum surface density is enough to have a functionally 
O ! diffcrcntiable Hamiltonian, thus providing definition of conserved quantities for the 

quasi-local regions. If on the boundary the evolution vector normal to the 2-surfacc is 



£h . chosen to be proportional to the dual expansion vector, we obtain a generalization of 

the Hawking energy associated with a generalized Kodama vector. This vector plays 
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the role for the stationary untrapped boundary conditions which the stationary Killing 
vector plays for stationary black holes. When the dual expansion vector is null, the 
' boundary conditions reduce to the ones given by the non-expanding horizons and the 

0^ ■ null trapping horizons. 
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1. Introduction 



Traditional description of black holes in terms of event horizons is inadequate for many 
physical applications, especially in the cases of non-stationary spacetimes. Quasi-local 
notions of trapped and marginal surfaces have now been found to be more useful 
for these cases within the framework of trapping, isolated, and dynamical horizons 
[H El El SI El El [TJ, [HI El OH EH]- These frameworks enable one to significantly extend 
the laws of black hole mechanics to the dynamical regime with the associated notions of 
energy, angular momentum and their fluxes, and have been applied to several problems 
in mathematical general relativity, numerical relativity, and quantum gravity |12j . These 
progresses on black hole dynamics lead to a question whether we can generalize the 
conservation laws for isolated and dynamical trapping horizons to general untrapped 
regions so that we can study the change of energy, angular momentum and their fluxes 
for untrapped strong gravitating systems, e.g., before the black hole horizon was formed. 

The question of how to define energy and angular momentum for untrapped surfaces 
has been raised for a while in searching for the "quasi-local energy-momentum and 
angular momentum" [13] . The goal has been to find a suitable definition of total energy- 
momentum and angular-momentum, surrounded by a spacelike 2-surface S, with § 2 
topology, in 4-dimensional spacetime M. The construction is quasi-local in the sense 
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that it refers only to the geometry of S (intrinsic metric, first fundamental form), the 
extrinsic curvatures (second fundamental forms) and the connection 1-forms on the 
normal bundle (normal fundamental forms) for its embedding in M. 

A systematic way to study conserved quantities is through the Hamiltonian. For 
the existence of a functionally differentiable Hamiltonian for General Relativity, it is 
necessary to impose suitable boundary conditions. The allowed boundary conditions for 
finite spatial 2-surfaces were studied previously with Dirichlet and Neumann boundary 
conditions [HJ [15]. This extends the requirement of the functional differentiability of 
the Hamiltonian, considered first by Regge and Teitelboim [16] for spatial infinity to the 
finite spatial 2-surfaces. Especially interesting boundary conditions are the Dirichlet 
boundary conditions which fix the induced metric on the 2-surfaces. 

The conditions were recently relaxed by Szabados who showed that instead of fixing 
the full induced metric on the boundary, fixing the area element is enough. Together with 
conditions that the lapse is vanishing on and the shift is tangent to the boundary and is 
divergence free, he showed that the Hamiltonian is functionally differentiable, and that 
in the large sphere limit the conserved quantities derived from the Hamiltonian behave 
as the spatial components of the total angular momentum [17]. This leads to a question 
whether or not the condition of vanishing lapse can be relaxed, and what additional 
conditions should be imposed so that a functionally differentiable Hamiltonian can 
provide a definition of energy-momentum for the quasi-local region. 

In this paper, we propose a class of quasi-local boundary conditions for canonical 
general relativity. It is shown that for untrapped or marginal surfaces, fixing the area 
element on the 2-surface S (rather than the induced 2-metric) is enough to have a 
functionally differentiable Hamiltonian, thus providing definition of conserved quantities 
for the quasi-local regions that allows the geometry outside to be dynamical and 
admit gravitational and other radiation. For gravitating systems including angular 
momentum, we further fix the " angular momentum surface density 11 at boundary to 
obtain a generalized definition of quasi-local energy including angular momentum. 
These boundary conditions characterize the equilibrium situation for regions bounded 
by untrapped or marginal surfaces. 

The evolution vector can be chosen freely. If on the boundary the evolution vector 
normal to the 2-surface is chosen to be proportional to the dual expansion vector, we 
obtain a generalization of the Hawking energy associated with a generalized Kodama 
vector. This vector plays the role for these stationary untrapped boundary conditions 
which the stationary Killing vector plays for stationary black holes. When the dual 
expansion vector is null, the boundary conditions reduce to the ones given by null 
trapping horizons and non-expanding horizons. 

We begin with the geometry of an untrapped 2-surface S embedded in a 4- 
dimensional spacetime M. Introduce a set of orthonormal vectors eo,ei,e2, e$ adapted 
to the 2-surface S, with eo and e\ being the set of timelike and spacelike unit normals 
to S and ca = [ea, being tangent to S. The extrinsic curvatures of S with 
respect to e and e± directions are given by k(e )AB = 5 , (e_B,V a^o) = —^oba and 
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k(ei)AB — g( e B,^A e i) — — Fiba- The connection 1-forms in the normal bundle are 
given by zua = ^(ei, V^eo) = — Tqia- Here T IJK = — g(ej,V k^i) are Ricci rotation 
coefficients. Under a boost transformation of eo and e±, 

e'o = eo cosh u + e\ sinh u, 

e'i = eo sinh u + t\ coshw, (1) 
they transform as 

fc(e'o) ab = k(e ) A B cosh u + k(ei) A B sinh u, 
k{e'i)AB = k(e )AB sinh u + k(ei)AB cosh u, 

hj'a = — Va«- (2) 

The expansion vector H, and the dual expansion vector H± are defined with the 
trace of the extrinsic curvatures k(e ) and k(ei), 

H = k(e 1 )e 1 - k(e )e , (3) 
H± = &(ei)e - k(e )e 1 , (4) 

where H and H± are also known as the mean curvature vector and the dual mean 
curvature vector respectively. These vectors are independent of choice of normal frames 
for the 2-surface. They are invariant under the boost transformation ([T|). Thus, they 
depend only on the 2-surface S and constitute a set of natural normal vectors for S [18J. 

A 2-surface S is trappedii k(ei) 2 > k(eo) 2 , untrappedii k(ei) 2 < k(eo) 2 , and marginal 
if fc(ei) 2 = k(e ) 2 , everywhere on S. Untrapped surfaces are also called mean convex 
surfaces. The dual expansion vector H± is always timelike for untrapped surfaces, null 
for marginal surfaces [HI [TH H3] . Note that, on S, the trace of the extrinsic curvature 
is zero along the direction of the dual expansion vector, i.e. k(H±)\s = 0, thus we have 

£h ± q\s = 0, (5) 
where g is the area element of S. 



2. Conserved quantities derived from a hamiltonian 

For a general diffeomorphism-invariant field theory in four dimensions with a Lagrangian 
4-form £(<p), where <fi denotes an arbitrary collection of dynamical fields, the field 
equations, £ = 0, are obtained by computing the first variation of the Lagrangian, 

5C = dQ((f),5<j)) + £6<t>, (6) 

where 0(0, £0) is the symplectic potential 3-form. For any diffeomorphism generated 
by a smooth vector field £, 

di^C = £^C = dQ((f>, £^<p) + ££^<p, (7) 

where £^ denotes the Lie derivative and i^ is the inner product, one can define a 
conserved Noether current 3-form J(£) by 

J(O = e(0,^0)-z e £(0), (8) 
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such that the Noether current J(£) is closed (dJ(£) = — ££^0 ~ 0) when the field 
equations are satisfied. Locally there exist a 2- form (called the Noether charge) 
such that J(£) = dQ(£). The variation of the Noether current 3-form is given by 

SJ(0 = u(<f>, 50, £(</>) + d(i^(<f>, 50)), (9) 

where u is the symplectic current 3-form defined by 

u((f>, 5 l( j), 62(f)) = 5^(0, 5 2 0) - 5 2 6(0, ^0). (10) 

Its integral over a 3-surface S defines the presymplectic form Q. If the presymplectic 
form is a total variation 

ft(0,50, £^) = J u(</>,6<f>,£t<i>) = 5H(£), (11) 

for some function H(£) on the field space, then H(£) is conserved along £, i.e. we 
have £^H(^) = 0. The function H(£) is called the Hamiltonian conjugate to £ 
[j~9l 1201 1211 1221 [23] . Note that on shell, the presymplectic form is given by 

/ u{<j>, 50, £ 5 0) = / 5J(£) - d{izQ) = I 6Q(£) - ^6. (12) 

For General Relativity we begin with the Hilbert action, 

S = J C = J R ab A *($a A (13) 

where i? ab = <fT ab + r a c A r cfe is the curvature 2-form constructed by the connection 
1-form T ab , *($ a A -tf 6 ) = \e ab cd -d c A and g> = r] ab § a <g> s § b is the metric, where 
— diag(— 1, 1, 1, 1) and d 0, is the orthonormal frame 1-form field. A variation of the 
Lagrangian gives, 

5C = R ab A e abcd $ c A 8$ d + 5T ab A D * (0 a A 6 ) 

+ d(5r l "A*(tf a A^)), (14) 

which identifies the symplectic potential 0(0,50) = 5T ab A *($ a A The Noether 
current 3-form is given by 

AO = dQ(0 = d[^T ab ) * (0 a A 4)], (15) 

where is the Noether charge 2-form and we assume the field equations are satisfied. 
The presymplectic form is given by 



/ u = I 8Q{0 -i t Q=j> d(0 + I C, 

t/S J S J S J s 



(0, (16) 



where, for convenience, we define 



and 



C 1 (0 = ^ ab 5(e abcd ^ c A^ d ), (17) 



C 2 (0=t^ c A5T ab Ae abcd $ d 
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Let us first expand C\ to the normal and tangent components, 

I d(£) = I isT 01 6(e 01AB $ A A i!) B ) + d) kT AB 5(e AB01 ^ A tf 1 ) 

J S J S J s 

+ 2 I i^ 0A 8(e AiB^ A d B ) + 2 / ^r 1A 5( eiA0B ^° A ^).(19) 

Let Ps be the projection onto S, we have Ps^ = Ps^ 1 = 0. Moreover, assuming 
<5(Ps) = 0, we have 

P S S$° = 5{P S $°) - {5P S )$° = 0, (20) 

and likewise PsS^ 1 = ^Ps?? 1 ) — (^Pj)?? 1 = 0. Thus the term with C\ reduces to 

/d(£) = 2 /^r 01 ^, (21) 

where = |(eoiAB$ A A i!) B ) is the area element on S. 

Similarly for C*2, by a projection onto the 2-surface S, it is straightforward to show 
that the following identity holds on S for any vector field V = V a e a : 

V c 5Y ab A e abcd ¥\ s = -2 g {V°5k( ei ) + VHk^) - V A 5w A )\ s . (22) 

Using this identity, we obtain 

/ c 2 {£) = -<[>2q (i^°Sk( ei ) + i^Skfo) - i^ A 5w A ) . (23) 
Js Js 

By (12T|) and (!23j) . the full symplectic form is, 



c^ = ^ C a (£) + C a (£) 

= / 2 z f r 01 5^ - / 2^ (i £ o <fjfe(ei) + z^tfA^eo) - i^5w A ) . (24) 
Js Js 

This is the key equation for our discussion. 
3. Energy 

We assume that the timelike (or null) vector £ is fixed on S, 

<^£|s — 0) (Boundary Condition I), (25) 
moreover, we assume 

$q\s — 0, (Boundary Condition II), (26) 

i.e. the area element of the 2-surface g is fixed. 

We first consider a special case by assuming that, on S, i^ A \s = (this condition 
will be relaxed in the next section). Using the boundary conditions (I, II), 

C 2 (£) = -S <f git (^k(e l )+^ 1 k(e )) , (27) 

S Js 

is a total variation, and 

i"d(£) = 2 li^Hg = 0. (28) 
Js Js 
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Thus with boundary conditions (I) and (II), the full symplectic form, 

/ u = -6 i (# fc(eO + ^k(eo)) = 6E H {£), (29) 
is a total variation. The Hamiltonian E#(£) associated with the vector £ is given by 

MO = £ (/(ff) - k (#°k(ei) + ^Heo))) Q. (30) 

where f(g) is a function of the area element. The evolution vector £ can be chosen 
freely. By choosing £ to be timelike or spacelike vector, the Hamiltonian gives energy 
or momentum respectively. Here we assume £ is timelike. 

Note that because of equation (jSJ), the boundary conditions (I) and (II) are both 
satisfied if we replace the variation with the Lie derivative with respect to £ (5 = 
and assume that on S, the evolution vector £ is given by 

Z\s = h{g)H ± , (31) 

where H± is the dual expansion vector and h{g) is a function of the area element on the 
(untrapped or marginal) 2-surface S. 

The boundary condition II implies that the area is conserved along the dual 
expansion vector. Since the area of each cross section does not change along the 
dual mean curvature vector direction, we called it the stationary untrapped boundary 
conditions. 

Note that because of equation fl3TT) . £|g = h(g)H±, so by equation (j4]), we have, on 
S, i^° = h(g)k(ei), i^fi 1 = —h(g)k(eo), and i^ A = 0. Thus the Hamiltonian E#(£) 
associated with the vector £ is given by 

Etf(£) = ^ (f(g) - h(g)H 2 ) g. (32) 

The free functions of the area element, f(g) and h(g), can be chosen such that the 
expression gives ADM mass at spatial infinity and irreducible mass at marginal surfaces 
H = 0. This can be done by letting f(g) to be l/(87rR) and let h(g) to be R/(327r), 
where K is the area radius given by 



1 

47T 



g. (33) 



s 



This leads to the energy expression, 



^^ft 1 -!^^)' < 34) 

which is precisely the Hawking energy [24J . 
4. Angular momentum 

The equation ff24l can also be used to define angular momentum. Let £|g = (87r) -1 ^ be 
a vector tangent to S satisfying i^°\s = and i^p^ls = 0, and 

5ip\s = 0, (Boundary Condition I'), (35) 
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and 

zua5$ a \s = 0, (Boundary Condition A), (36) 



then 



u = —5 & i) A w A Q = W), (37) 
8n 



is a total variation, where ip = i^fi . The Hamiltonian associated with ip is the angular 
momentum given by 

J(V>) = ?T ( f3Q=ir ( f Q, (38) 

8tt J s 8tt J 5 

with the "angular momentum surface density" given by j = tp A WA- 

Because of the gauge freedom in choosing the normal fundamental forms w A 
(equation ([2j)), the definition of angular momentum is not unique. A further condition 
which makes it unique is 

$tpQ\s = £ipQ\s = 0, (Boundary Condition II'). (39) 

Here ip generates a symmetry of the area form rather than the whole metric. This implies 
that ip has vanishing transverse divergence V A ip A = 0. Under the gauge freedom for 
the normal fundamental forms wa | — > — V^m, the angular momentum formula 

& ip A zu A i-> I i) A {vJ A - V a m) = i 4> a tu a - V A {^ A u) = I i\) A w A , (40) 
Js Js Js Js 

is invariant. Thus the angular momentum formula is uniquely defined on S (up to the 

choice of ip A ). This condition is also given in [6j [TOj [11] for dynamical trapping 

horizons and was also used by Szabados [17j as a Hamiltonian boundary condition for 

quasi-local angular momentum. 

We can now extend the stationary untrapped boundary conditions to include 

angular momentum. Similar to the discussion in the previous section, we require that 

$£\s — (Boundary Condition I) and Sg\s = (Boundary Condition II) are satisfied. 

In addition, we further require that Boundary Condition A is satisfied. This implies 

that the angular momentum surface density j is fixed on S, i.e. 

$j\s — 0, (Boundary Condition III). (41) 

The full symplectic form is then given by applying the boundary conditions to the 
equation (T5^ 



w(0,<ty,£ e 0)=5E(£), (42) 
where the Hamiltonian associated with £ is now given by 

E(0 = j (f(Q,j) ~ k (0°*(ei) + ^(eo))) Q- (43) 
A natural choice of the evolution vector on S is 

Z\s = h(Q,j)H ± -£l{Qj)iP, (44) 
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which is assumed to be timelike or null. The free functions h(g,j) and Q(g,j) (angular 
speed) are now functions of the area element q and the angular momentum surface 
density j. 

Now we replace the variation S by the Lie derivative £g, the boundary condition 
(I) is automatically satisfied (£^£ = 0). By equation ( 1441) . the boundary condition (II) 
is satisfied with the following condition on S: 

^Q\s = 0. (45) 

The boundary condition (III) is satisfied if the vector field ip is transported by H±, 

£^H ± \ S = -£h ± ?P\s = 0. (46) 

The conditions fj4"5]) and (146]) are consistent with the ones in [6l [9j [TOj, [11] for dynamical 
trapping horizons. These conditions then imply that 

£f£|s = 0, £ 5 £?|s = 0, £ € j|s = 0, (47) 

the stationary untrapped boundary conditions (I, II, III) are satisfied. 

By equation (l4~3l) . the Hamiltonian associated with the evolution vector (144p is then 
given by 

E(0= <f (f(Q,j)-h(Q,j)H 2 ))g. (48) 

By requiring that the energy expression gives the standard value for Kerr black holes 
at marginal surface H = 0, 

= ^3!, (49) 

we obtain 



The other free function h(g,j) can be chosen such that the energy is proportional to the 
Hawking energy and gives ADM mass at spatial infinity, this implies that 

VR 4 + 4JP 

h{Q > 3) = ^2^-- (51) 
Then the Hamiltonian associated with £, 



= ^EpE f! - 1 f „> Q . , (52) 



2M V 167r 
provides a suitable choice for the energy expression. 
Note that the vector on S\ 



{ | S = v™ Hl _ W) , (53) 
327rlR 

reduces to the Kodama vector 

R 

^Kodama = (54) 
OZ7T 

in spherically symmetric spacetimes [25]. Thus equation (|53l) is a generalized Kodama 
vector for non-spherically symmetric spacetimes. 
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5. Discussion 

In summary, for quasi-local regions bounded by an untrapped or marginal 2-surface, 
a functional differentiable Hamiltonian can be denned with the boundary conditions 
which fix the evolution vector (equation (25)), the area element (equation (26)) and the 
angular momentum surface density (equation (41)) on the 2-surface. As a consequence, 
a preferred expression of quasi-local energy for these "stationary untrapped boundary 
conditions" is given, which allows the geometry outside to be dynamical and to admit 
gravitational and other radiation. The quasi-local energy expression (equation (52)) 
generalizes the Hawking energy to include angular momentum. On the boundary, the 
evolution vector associated with this expression is a generalization of the Kodama vector 
(equation (53)). These results generalize the conserved quantities for Isolated Horizons 
so as to provide covariant conserved quantities for general untrapped regions. 

In equation (31), the evolution vector on the boundary is chosen to be proportional 
to the dual expansion vector. Alternatively, we can use the unit dual expansion vector. 
If H is not null {\H\ ^ 0), there is a set of uniquely determined unit normal vectors for 
the 2-surface, given by 



where \H\ — \JH ■ H = ^/k(ei) 2 — k(eo) 2 . If we choose the evolution vector £ such that 
£|s — eo, then by equation ( 1301) . this leads to the energy expression 



Here we have only one free function f(g). A natural requirement is that the expression 
should give ADM energy at spatial infinity. Assuming that we can embed the 2-surface 
isometrically into Minkowski spacetime, let fco(ei) be the trace of extrinsic curvature with 
respect to e±, for the 2-surface in Minkowski spacetime, then the choice f(g) = k (ei) 
gives the Kijowski-Liu-Yau quasi-local energy [261 [27]. In the special fc(eo) = slice, the 
Kijowski-Liu-Yau energy equals the Brown- York energy [281 EH] - There are nice positive 
energy theorems for these expressions [27, 30] . However, eo fails to be defined in the null 
case. Moreover, the requirement of isometric embedding assumes that the full induced 
2-metric is fixed (rather than just the area element). Thus it seems that this choice is 
not suitable for the cases involving dynamical black holes. 

The stationary untrapped boundary conditions can be compared with the Isolated 
Horizon boundary conditions. The untrapped surface S together with the associated 
evolution vector on the boundary £|g constitute a timelike hypersurface. In the 
limit when the dual expansion vector H± is null, S reduces to a marginal surface, 
the hypersurface reduces to a non-expanding horizon (null 3-surfaces in spacetime on 
which the expansion with respect to any null normal is zero) [12]. This suggests a 
generalization of the non-expanding horizons to "stationary untrapped hypersurfaces" 
(with the boundary conditions I, II, III) describing the equilibrium states of untrapped 




(56) 
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surfaces with conservation of the area (with area radius R), the energy E, and the 
angular momentum J. 
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